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, Astrophysical jets are associated with the formation of young stars of 

qq * all masses, stellar and massive black holes, and perhaps even with the 

formation of massive planets. Their role in the formation of planets, 
stars, and galaxies is increasingly appreciated and probably reflects a 
deep connection between the accretion flows - by which stars and black 
holes may be formed - and the efficiency by which magnetic torques can 
remove angular momentum from such flows. We compare the properties 
and physics of jets in both non-relativistic and relativistic systems and 
trace, by means of theoretical argument and numerical simulations, the 
physical connections between these different phenomena. We discuss 
the properties of jets from young stars and black holes, give some basic 
theoretical results that underpin the origin of jets in these systems, and 
then show results of recent simulations on jet production in collapsing 
star-forming cores as well as from jets around rotating Kerr black holes. 
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1.1 Introduction 



The goal of this book, to explore structure formation in the cosmos and 
the physical linkage of astrophysical phenomena on different physical 
scales, is both timely and important. The emergence of multi- wavelength 
astronomy in the late 20th century with its unprecedented ground and 
space-based observatories, as well as the arrival of powerful new com- 
putational capabilities and numerical codes, have opened up unantici- 
pated new vistas in understanding how planets, stars, and galaxies form. 
Galaxy formation has turned out to be a complex problem which requires 
a deep understanding of star formation on galactic scales and how it feeds 
back on galactic gas dynamics. The discovery of significant numbers of 
massive black holes in galactic nuclei now suggests that most galaxies 
harbour million to billion mass holes. The formation of these monsters 
and the effects that they can have on the early evolution of galaxies is 
still not well understood however. Star formation studies have made 
huge inroads, but must still take a large step before we can truly under- 
stand the origin of stellar masses and star formation rates in molecular 
clouds and galaxies. The newly discovered planetary systems bear lit- 
tle relation to text-book models of how our solar system is believed to 
have formed. It is also abundantly clear that planet and star formation 
are intimately coupled through the physics of protostellar disks that are 
their common cradles. 

Astrophysical jets play an important role in the formation and evolu- 
tion of stars, black holes, and perhaps massive planets. They appear to 
be an inescapable multi-scale phenomenon that arises during structure 
formation. The enormous kinetic luminosity of many jets (often com- 
parable to the bolometric luminosity of the central sources) also implies 
that they could have important feedback effects on structure formation; 
from the scale of the giant lobes of radio galaxies stretched out across 
many Mpc scales in the IGM, to jets from stellar mass black holes, down 
to the stirring up of molecular gas on sub-pc to pc scales in regions of 
low mass star formation. Jets are observed to be ubiquitous during the 
process of star formation and associated with a very broad range of ob- 
jects; from brown dwarfs to B and perhaps even O stars. It is now well 
established that stellar mass black holes are associated with jets and that 
the properties of such " microquas ars" scale very well only with the mass 
of the black hole (e.g. see review, Mirabell . 2005 ). and scale naturally to 



the limit of massive black holes. In this regard, galactic jets are being 
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increasingly asso ciated with accreting massive black holes (e.g. review, 



Blandfordll200lh . 



In most cases, there is clear evidence that accretion disks are an es- 
sential part of the engine. The fact that outflows and jets are observed 
around nearly all astrophysical objects during the early stages of their ex- 
istence - a time where the central objects (young stars, massive planets, 
or black holes) undergo significant accretion from surrounding collapsing 
gas structures and/or associated disks - argues for a deep link between 
accretion and outflow. A large body of observations and theoretical mod- 
els increasingly suggests that such jets in all of these systems may be 
powered by the same mechanism, namel y, hydromagnetic winds driven 
off of magnetised ac c retion disks (see e.g. Liviol . 1999 : Konigl fc Pudrit 



2000; Pudritz et al. . 2007 . for reviews on this subj ect). The basic the- 



ory fo r hydromagnetic disk winds was worked out by Blandford fc Pavnel 
( 1982h in the context of models for AGN jets, but was quickly found to 
be very important for understanding jets in protostellar systems (PN83). 
For accreting systems, the gravitational binding energy that is released 
as gas accretes through a disk and onto a central object is the ultimate 
source of the energy for jets. Therefore, if this energy can be efficiently 
tapped and carried by the jet, then jet energies and other properties 
should simply scale with the depth of the gravitational potential well 
created by the central mass. The torques that can be exerted by out- 
flows upon the underlying disks were shown t o be much larger than can 
be produced by even stron gly turbulent disks (jPudritz fc Normanl ll983; 
Pelletier fc Pudrit zl . fl99l) . Thus, jets can also carry off most of the an- 



gular momentum underlying disks, thereby assisting in the growth of 
central objects - be they stars, black holes, or planets. (It should be 
noted that large scale spiral waves in accretion disks can also efficiently 
transport angular momentum through disks). Finally, the feedback of 
these powerful jets upon their surroundings can be important. In star 
forming systems, outflows may help to drive turbulence within cluster 
forming regions in molecular clouds, which would help to sustain such a 
region against global collapse. 

It must be noted that isolated magnetised, spinning bodies such as 
magnetised A stars can also drive drive outflows. The extraordinary 
Chanra X-ray observations of Crab pulsar show that it driving off a 
highly collimated, relativistic jet. The models that best describe the 
physics of all of these jets demonstrate that outflows can be driven and 
collimated from spinning bodies (be they disks, stars, or compact ob- 
jects). In protostellar jets, the slow spin of stars has been attributed 
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either to star-disk coupling (jKonigll . Il99lh or X- winds |Shu et alllioooh . 
but could also be e xplained as an accretio n powered stellar wind from 
the central star (e.g. jMatt fc Pudrita . 120051 ). For black hole systems, the 
highly relativistic component of jets may be associated with electrody- 

namic processes and currents that arise from spinning hol es in a magne- 

tised environment provided by the surrounding disk (e.g-. lBlandford fc Znajek . 



1977). 



It is impossible to do justice to the enormous body of excellent work 
on the physics of jets in these different systems. In this chapter, we give 
specific examples of these basic themes from our own research in these 
fields in the context of some of the basic literature. In §2 we review 
aspects of unified models for non-relativistic, and then relativistic jets 
in these systems. We follow this in §3 with a discussion of disk winds 
as an excellent candidate for a unified model, analysed particularly for 
protostellar jets. We follow this in §4 with simulations of jets during 
gravitational collapse and disk formation. We examine the role of feed- 
back of jets in the specific context of star formation in §5. We then move 
on to the basics of relativistic jets from compact objects as well as black 
holes (§6), and then simulations (§7) of MHD jets in the general rela- 
tivistic limit (GRMHD). We conclude with a comparison of the physics 
of relativistic vs non-relativistic jets. 



1.2 Jets in diverse systems 

Protostellar objects : Stars span about 4 decades in mass. Until 
recently, the study of outflows was restricted to sources that were easily 
detectable in millimetre surveys of molecular clouds. These studies show 
that outflows, over a vast range of stellar luminosities, have scalable 
physical properties. In particular the observed ratio of the momentum 
transport rate (or thrust) carried by the CO molecular outflow to the 
thrust th at can be provided by the bolometric luminosity of the central 
star (e.g., Cabrit fc Bertoutl 1992f k scales as 

^outflow/^rad = 250( J L bo i/10 3 i Q )-°- 3 , (1.1) 

This relation has been confirmed and extended by the analysis of data 
from over 3 90 outflow s , rang ing over six decades up to 10 6 Lq in stellar 



luminosity (|Wu et all l2004f ). This remarkable result shows that that 
jets from both low and high mass systems, in spite of the enormous 
differences in the radiation fields from the central stars, are probably 
driven by a single, non-radiative, mechanism. 
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Jets are observed to have a variety of structures and time-dependent 
behaviour - from internal shocks and moving knots to syst ems of bow 



shock s that suggest long-time episodic outbursts (e.g., review iBallv et al 



2007). They show wiggles and often have cork-screw like structure, 



suggesting the presence either of jet precession, the operation of non- 
axisymmetric kink modes, or both. Numerical simulations have recently 
shown that in spite of the fact that jets are predicted to have a dominant 

toroidal magne tic field, they are nevertheless stable against such nonlin- 

ear m odes (e.g. JOuved et al.l . l2003tlNakamura fc Meierl . 12004 : iKigure fc Shibata 
2005). The nonlinear saturation of the unstable modes is achieved by 



the natural regulation of the jet velocity to values near the local Alfven 
speed. A second point is that jets will have some poloidal field along 
their "backbone" and this too prevents a jet from falling apart. 

Two types of theories that have been proposed to explain protostel- 
lar jets, both of w hich are hydromagnetic; disk winds (e.g., review 
Pudritz et~all l2007h . or the X-wind (e.g., review Ishang et all . l2007t ). 



The latter model posits the interaction between a spinning magnetised 
stellar magnetosphere and the inner edge of an accretion disk as the 
origin of jets, while the former envisages jets as arising from large parts 
of magnetised disks. One of the best ways of testing such theories is 
by measuring jet rotation - which is a measure of how much angular 
momentum that a jet can extract from its source. Significant observa- 
tional progress has been made on this problem lately with the discovery 
of the rotation of jets from T-Tauri stars. The angular momentum that 
is observed to be carried by these rotating flows (e.g. DG Tau) is a con- 
sider able fraction of t he excess disk angular momentum - from 60-100% 
(e.g.. iBacciottil . 120041 ). which is consistent with the high extraction effi- 
ciency that is predicted by the theoretical models. Another important 
result is that there is too much angular momentum in the observed jet 
to be accounted for the spinning star or very innermost region of the 
disk as predicted by the X-wind model. The result is well explained by 
the disk wind model which predicts that the angular momentum derives 
from large reaches of th e disk, typically from a region as large as 0.5 AU 



( Anderson et all l2003h 



Jovian planets : The striking similarity of the system of Galilean 
moons around Jupiter, with the sequence of planets in our solar sys- 
tem, has long suggested that these moons ma y have formed through a 
subdisk around Jupiter (jMohantv et al.l . l2007h . Recent hydrodynamical 
simulations of circumstellar accretion disks containing and building up 
an orbiting protoplanetary core have numerically proven the existence 
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of a cir cum-planetary su b-disk in almost Keplerian rotation close to the 



planet (jKlev et all 120011 ) 



M cp = 6x 10- 5 M Jup yv 



. The accretion rate of these sub-disks is about 
_1 and is confirmed by many independent sim- 



ulations. With that, the circum-planetary disk temperature may reach 
values up to 2000 K indicating a sufficient degree of ionisation for matter- 
field coupling and would also allow for strong equipartition field strength 
|Fendtl . [200a) . i should be possible, therefore, for lower luminosity jets 
to be launched from the disks around Jovian mass planets. 

The possibility of a planetary scale MHD outflow, similar to the larger 
scale YSO disk winds, is indeed quite likely because: (i) the numerically 
established existence of circum-planetary disks is a natural feature of 
the formation of massive planets; and (ii) the feasibi lity of a large-scale 
magnetic field in the protoplanetary environment ( Quillen fc Tri lling. 
19981 : Fendt . 2003 ). One may show, moreover, that the outflow velocity 



is of the order of the escape speed fo r the protoplanet, at about 60 km s 



|Fendtl . liooi iMachida et allbood ). 



Black holes : Relativistic jets have been observed or postulated 
in va rious astrophysical objects, including a ctive galactic nuclei (AGNs) 
(e.g., 



Urrv fc Padovani[ 1995 : Ferrari . 1998 ). microquasars in our galaxy 
Mirabel fc Rodriguez! 1998 ). and gamma-ray bursts (GRBs) (e.g., 



Piranl 120051 ). Table [LT1 summarises the features or relativistic jets in 
different large scale sources and demonstrate the wide range in power 
and Lorentz factors achieved by these systems. 



Table 1.1. The sources and features of large scales /relativistic jets 



Source 




L 






r 


AGNs / Quasars 


10 43 


10 48 


erg 


s- 1 


~ 10 


/i-quasars 


10 3S - 


■10 40 


erg 


s" 1 


1-10 


GRBs 


10 51 - 


10 52 


erg 


s" 1 


10 2 -10 4 



In the commonly accept ed standard model of large scale/relativistic 



jets (jBegelman et allll984l ). flow velocities as large as 99% of the speed 



of light (in some cases even beyond) are required to explain the appar- 
ent superluminal motion observed in many of these sources (see Table 
ll.ljl . Later, considerations of stationary MH D flows have revealed that 
relativistic jet s must be strongly magnetised (Michell . 19691 : Camenzindl . 
1986tlLi et al] . ll992HFendt fc Camenzindl . Il996h . In that case, the avail- 
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able magnetic energy can be transferee! over a small amount of mass with 
high kinetic energy. 

Here too, the most promising mechanisms for producing the relativis- 
tic jets involve magnetohydrodynamic centrifugal acceleration and/or 
magnetic pressure driven acceleration from an acc retion disk around 
the compact objects (e.g.. lBlandford fc PayneL Il982j ). or involve the ex- 
traction of rotating energy from a rotating black hole (jPenrosel [1969; 



Blandford k, Znaiekl . I1977N . These models have been ap plied to explain 



jets fe atures in the galactic micro-quasars GRS 1915+105 ( Mirabel fc Rodriguez! 



1994) and GRO J1655-40 (jTingav et al.l . ll995l ). or a rotating super mas- 



sive black hole in an active galactic nucleus, which is fed by interstellar 
gas and gas from tidally disrupted stars. In general these studies re- 
quire solving special relativistic MHD (SRMHD) or general relativistic 
MHD (GRMHD) equations and often require sophisticated numerical 
codes. We describe the basic theory behind relativistic MHD jets and 
summarise recent GRMHD simulations of jets emanating from the vicin- 
ity of accreting black holes (see §6 and 7) after we examine the general 
theory of outflows from magnetised spinning disks and objects.. 



1.3 Theory of disk winds 

Given the highly nonlinear behaviour of the force balance equation for 
jets (the so-called Grad-Shafranov equation), theoretical work has fo- 
cused on tractable and idealised time-independent, and axisymmetric or 
self-similar models (e.g., BP82) of various kinds. We briefly summa rise 



the theory, below (see details in lPudritz et al.l . l2007t iPudrita . 120041) . 



1.3.0.1 Conservation Laws and Jet Kinematics 
Conservation laws play a significant role in understanding astrophysical 
jets. This is because whatever the details, conservation laws strongly 
constrain the flux of mass, angular momentum, and energy. What can- 
not be constrained by these laws will depend on the general physics of 
the disks such as on how matter is loaded onto field lines. 

Jet dynamics can be described by the time-dependent, equations of 
ideal MHD. The evolution of a magnetised, rotating system that is 
threaded by a large-scale field B involves (i) the continuity equation 
for a conducting gas of density p moving at velocity v (which includes 
turbulence); (ii) the equation of motion for the gas which undergoes 
pressure (p), gravitational (with potential $), and Lorentz forces; (hi) 
the induction equation for the evolution of the magnetic field in the 
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moving gas where the current density is j = (c/47r)V x B; (iv) the en- 
ergy equation, where e is the internal energy per unit mass; and, (v) the 
absence ol magnetic monopoles. These are written as: 



dp 

dt 



V.(pv) = (1.2) 



^ + (v.V)v ] + Vp + pV$ - = (1.3) 

at J c 

— -Vx(vxB) = (1.4) 

P fir + (v-V)e^) +p(V.v) = (1.5) 



dt 



V.B = (1.6) 



Progress can be made by restricting the analysis to stationary, as well 
as 2D (axisymmetric) flows, from which the conservation laws follow. 
It is useful to decompose vector quantities into poloidal and toroidal 
components (e.g. magnetic field B = B p + B^e^). In axisymmetric con- 
ditions, the poloidal field B p can be derived from a single scalar potential 
a{r, z) whose individual values, a = const, define the surfaces of constant 
magnetic flux in the outflow and ca n be specified at the surface of the 
disk (e.g JPelletier fc Pudritzl . ll992l PP92). 



Conservation of mass and magnetic flux along a field line can 
be combined into a single function k that is called the "mass load" of 
the wind which is a constant along a magnetic field line; 

pvp = fcB p . (1.7) 

This function represents the mass load per unit time, per unit magnetic 
flux of the wind. For axisymmetric flows, its value is preserved on each 
ring of field lines emanating from the accretion disk. Its value on each 
field line is determined by physical conditions - including dissipative 
processes - near the disk surface. It may be more revealingly recast as 

fc(a) = b; = -a* - ' (L8) 

where dM w is the mass flow rate through an annulus of cross-sectional 
area dA through the wind and d^ is the amount of poloidal magnetic 
flux threading through this same annulus. The mass load profile is a 
function of the footpoint radius ro of the wind on the disk. 

The toroidal field in rotating flows derives from the induction equa- 
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tion; 

B 4 , = £(v4,-n r), (1.9) 

where fi is the angular velocity of the disk at the mid-plane. This 
result shows that toroidal fields in the jet are formed by winding up the 
field from the source. Their strength also depends on the mass loading 
as well as the jet density. Denser winds should have stronger toroidal 
fields. We note however, that the density does itself depend on the 
value of k. Equation (9) also suggests that at higher mass loads, one has 
lower toroidal field strengths. This can be reconciled however, since it 
can be shown from the conservation laws (see below) that the value of 
k is related to the density of the outflow at the Alfven point on a field 
line; k — (pa/^k) 1 / 2 (e.g., PP92). Thus, higher mass loads correspond 
to denser winds and when this is substituted into equation (9), we see 
that this also implies stronger toroidal fields. 

Conservation of angular momentum along each field line leads 
to the conserved angular momentum per unit mass; 

tBa, 

1(a) = r«0 - — j- = const. (1-10) 

The form for I reveals that the total angular momentum is carried by 
both the rotating gas (first term) as well by the twisted field (second 
term), the relative proportion being determined by the mass load. 

The value of 1(a) that is transported along each field line is fixed by 
the position of the Alfven point in the flow, where the poloidal flow 
speed reaches the Alfven speed for the first time (wa = 1)- It is easy to 
show that the value of the specific angular momentum is; 

1(a) = Q r 2 A = (r A /r Q ) 2 l . (1.11) 

where Iq = VKflfo — ^o?"o i s the specific angular momentum of a Kep- 
lerian disk. For a field line starting at a point rg on the rotor (disk in 
our case), the Alfven radius is r A (r ) and constitutes a lever arm for the 
flow. The result shows that the angular momentum per unit mass that 
is being extracted from the disk by the outflow is a factor of (r A /r ) 2 
greater than it is for gas in the disk. For typical lever arms, one particle 
in the outflow can carry the angular momentum of ten of its fellows left 
behind in the disk. 

Conservation of energy along a field line is expressed as a gener- 
alised version of Bernoulli's theorem (this may be derived by taking the 
dot product of the equation of motion with B p ). Thus, there is a specific 
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energy e(a) that is a constant along field lines, which may be found in 
many papers (e.g., BP82 and PP92). Since the terminal speed v p = 
of the disk wind is much greater than its rotational speed, and for cold 
flows, the pressure may also be ignored, one finds the result: 

Voo ^ 2 1/2 Sl r A = (r A /r )v csc . Q . (1-12) 

There are three important consequences for jet kinematics here; (i) 
that the terminal speed exceeds the local escape speed from its launch 
point on the disk by the lever arm ratio; (ii) the terminal speed scales 
with the Kepler speed as a function of radius, so that the flow will have 
an onion-like layering of velocities, the largest inside, and the smallest 
on the larger scales, as seen in the observations; and (iii) that the ter- 
minal speed depends on the depth of the local gravitational well at the 
footpoint of the flow - implying that it is essentially scalable to flows 
from disks around YSOs of any mass and therefore universal. 

Another useful form of the conservation laws is the combination of 
energy and angular momentum conservation to produce a new constant 
along a field li ne (e.g., PP92) ; 7(a) = e(a) — flol(a). This expression 
has been used ( Anderson et al. . 2003) to deduce the rotation rate of the 



launch region on the Kepler disk, where the observed jet rotation speed 
is v 4>,oo at a radius r^ and which is moving in the poloidal direction with 
a jet speed of v p .oo. Evaluating j for a cold jet at infinity and noting 
that its value (calculated at the foot point) is j(ao) = — (3/2)w^- , one 
solves for the Kepler rotation at the point on the disk where this flow 
was launched: 

^Vp iOC / {2v 4>yO0 r 0O ) . (1.13) 
When applied to the obser ved rotation of the Lar ge Velocity Component 



(LVC) of the jet DG Tau (jBacciotti et all 120021 ). this yields a range of 



disk radii for the observed rotating material in the range of disk radii, 
0.3-4 AU, and the magnetic lever arm is r^/r^ — 1.8-2.6. 

1.3.0.2 Angular Momentum Extraction 

How much angular momentum can such a wind extract from the disk? 
The angular momentum equation for the accretion disk undergoing an 
external magnetic torque may be written: 



(1.14) 
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where we have ignored transport by MM turbulence or spiral waves. By 
using the relation between poloidal field and outflow on the one hand, 
as well as the link between the toroidal field and rotation of the disk 
on the other, the angular momentum equation for the disk yields one of 
the most profound scaling relations in disk wind theory - namely - the 
link between disk accretion and mass outflow rate (see KPOO, PP92 for 
details): 

M a ~ (r A /r ) 2 M w . (1.15) 

The observationally well known result that in many systems, M w /M a ~ 
0.1 is a consequence of the fact that lever arms are often found in nu- 
merical and theoretical work to be ta/tq ~ 3 - the observations of DG 
Tau being a perfect example. 

1.3.0.3 Jet Power and Universality 

These results can be directly connected to the observations of momen- 
tum and energy transport in the molecular outflows. Consider the total 
mech anical power t hat is carried by the jet, which may be written as 



(e.g.. |Pudritd . l2003h : 




L jet = \ £ dM w vl ~ ^1-^)^ li.cc. (1.16) 

This explains the observations of Class outflows wherein L w /Lboi — 
1/2, since the main luminosity of the central source at this time is due 
to accretion and not nuclear reactions. (The factor of 1/2 arises from 
the dissipation of some accretion energy as heat at the inner bound- 
ary). The ratio of wind to stellar luminosity decreases at later stages 
because the accretion luminosity becomes relatively small compared to 
the bolometric luminosity of the star as it nears the ZAMS. 

This result states that the wind luminosity taps the gravitational en- 
ergy release through accretion in the gravitational potential of the cen- 
tral object - and is a direct consequence of Bernoulli's theorem. This, 
and the previous results, imply that jets may be produced in any accret- 
ing system. The lowest mass outflow that has yet been observed corre- 
sponds to a proto-brown dwarf of luminosity ~ 0.09Lq, a stellar mass of 



only 20 — 45Mj up , and a very low mass disk < 10 4 M© (jBourke et al 
20051) . 



On very general grounds, disk winds ar e also likely to be active dur- 
ing massive star formation (e.g., Konigj 1999h . Such outflows may 



already start during the early collapse phase when the central YSO 
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still has only a fract i on of a solar mass (e.g.. iBaneriee fc Pudritzl . [2006; 
Baneriee fc Pudrit3 . 2007 ). Such early outflows may actually enhance 
the formation of massive stars via disk accretion by punching a hole in 
the infalling envelop a nd releasing the building radiation pressure (e.g., 



Krumholz ct al 



2005). 



1.3.0.4 Jet Collimation 
In the standard picture of hydromagnetic winds, collimation of an out- 
flow occurs because of the increasing toroidal magnetic field in the flow 
resulting from the inertia of the gas. Beyond the Alfven surface, equa- 
tion (8) shows that the ratio of the toroidal field to the poloidal field 
in the jet is of the order B^/Bp ~ r/r^ ^> 1, so that the field be- 
comes highly toroidal. In this situation, collimation is achieved by the 
tension force associated with the toroidal field which leads to a radi- 
ally inwards directed component of the Lorentz force (or "z-pinch"); 
-pLorentz.r — jzBif>- The stability of such systems is examined in the next 
section. 



In iHevvaerts fc Normanl (|1989f ) it was shown that two types of solu- 
tion are possible depending upon the asymptotic behaviour of the total 
current intensity in the jet; 



/ = 2tt 



j z (r',z')dr' = (c/2)rB . 



(1.17) 



In the limit that / — » as r — > 00, the field lines are paraboloids which 
fill space. On the other hand, if the current is finite in this limit, then the 
flow is collimated to cylinders. The collimation of a jet therefore depends 
upon its current distribution - and hence on the radial distribution of 
its toroidal field - which, as we saw earlier, depends on the mass load. 
Mass loading therefore must play a very important role in controlling 
jet collimation. 



It can be shown (jPudritz et all 12006. PRO) from this that jets should 
show different degrees of collimation depending on how they are mass 
loaded. As an example, neither the highly centrally concentrated, mag- 
netic field lines associated with t he initial split-mqnopole magnetic con- 
figuration used in simulations bv lRomanova et alJ (I1997|K nor the simi - 



lar field structure invoked in the X-wind (see review bv lShu et all 12000) 



should become collimated in this picture. On the other hand, less cen- 
trally (radially) concen trated magnetic co nfigurations such as the po- 
tential configuration of lOuved et all (|1997f ) and BP82 should collimate 
to cylinders. 
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Fig. 1.1. Large scale outflow (left panel, scale of hundreds of AU), and and 
small scale disk wind and jet formed (right panel, scale of a fraction of an AU) 
during the gravitational collapse of a magnetised B-E, rotating cloud core. 
Cross-sections through the disk and outflows are shown - the blue contour 
marks the Alfven surface. Snapshots taken of an Adapti ve Mesh calculation 
at abo ut 70,000 years into the collapse. [Adapted from IBaneriee fc Pudritd 
(2003)]. 



This result also explains the range of collimation that is observed for 
molecular outflows. Models for observed outflows fall into two general 
categories: the jet-driven bow shock picture, and a wind-driven shell 
picture in which the molecular gas is driven by an underlyi ng wide-angle 



wind c omponent such as given by the X-win d (see revi e w by Cabrit et al . 
( 1997 ). A survey of molecular outflows by Lee et al. ( 2000h found that 
both mechanisms are needed in order to explain the full set of systems 
observed. 



1.4 Gravitational collapse, disks, and the origin of outflows 

Jets are expected to be associated with gravitational collapse because 
disks are the result of the collapse of rotating molecular cloud cores. 
One of the first simula tions to sh o w how jets arise during gravitational 
collapse is the work of iTomisakal (j 19981 . 120021 ) . These authors used as 
initial conditions the collapse of a magnetised rotating filament (cylin- 
der) of gas and showed that this gave rise to the formation of a disk 
from which a centrifugally driven disk wind was produced. 

Banerjee & Pudritz re visited the problem o f colla psing magnetised 
molecular cloud cores in Banerjee fc PudritJ ( 20061 ) (low mass cores, 
BP06) and IBaneriee k Pudritzl ( 2007 ) (high mass cores, BP07). Here 
the initial cores are modelled as supercritical Bonnor-Ebert (B-E) spheres 



14 The role of jets in the formation of planets, stars, and galaxies 





Fig. 1.2. Magnetic field line structure, outflow and disk. The two 3D images 
show the magnetic field lines, isosurfaces of the outflow velocities and isosur- 
faces of the disk structure at the end of our simulation (t ~ 7 x 10 4 years) 



at the two different scales as shown in Fig. 11.11 The isosurfaces of the upper 
panel refer to velocities 0.18kms _1 (light red) and 0.34 km s" 1 (red) and a 
density of 2 x 1CT 16 gem" 3 
faces with velocities 0.6 km s 



jray) whereas the lower panel shows the isosur- 
1 (light red) a nd 2kms~ 1 (red) and the d ensity 



at 5.4 x 10" 9 gcm" d (gray). [Adapted from iBaneriee fc Pudritzl ([2006)]. 



with a slight spin. B-E spheres are well controlled initial setups for 
numerical experiments and th ere are plenty of observed cores which 



can fit by a B-E profile (see lLada et all 120071 . on a compilation of 



hydrostatic cores). Recent efforts have focused on understanding the 
evolution of magneti se d B-E spheres (e.g. Matsumoto fc Hanawai l2003t 
Machida et~aH . l2005l l?l: lHennebelle fc Tevssieiil2007tlHennebelle fc Fromane 



20081) . Whereas p urely hydrodynamic collapses of such ob j ects n ever 
show outflows (e.g., Foster fc Chevalier . 19931 : Baneriee et al. . 2004 ). the 
addition of a magnetic field produces them. 

Outflows are also implicated in playing a fundamental role in con- 
trolling the overall amount of mass that is incorporated into a new star 
from the original core. Observations of the so-called core mass func- 
tion (CMF) that describes the spectrum of core masses, show that it 
has nearly the identical shape as the initial mass function (IMF) that 
describes the spectrum of initial stellar masses - with one caveat. The 
CM F needs to be mulip lied by a factor of 1/3 (in one well studied case 



sec 



Alves et all 120071 ) in order to align them. The suggestion is that 



protostellar outflows are responsible for removing as much as 2/3 of a 
core's mass in the process of collapse and star formation. 

In BP06, the FLASH AMR code ijFryxell et all boooh was used to 
follow the gravitational collapse of magnetised molecular cloud cores. 
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Fig. 1.3. Snapshots of the central region of a collapsing magnetised massive 
cloud core (JWcore ~ 170 M©). The left panel shows the situation at t = 
1.45 x 10 4 years (1.08%) into the collapse and before the flow reversal and 
right panel shows the c onfiguration 188 years later when the outflow is clearly 
visible. [Adapted from iBaneriee fc Pudritd i|2007h ]. 



This code allowed the Jeans length to be resolved with at least 8 grid 
points throughout the colla pse calculations . Thi s ensures compliance 
with the Truelove criterion (jTruelove et al. . 1997 ) to prevent artificial 
fragmentation. 

The results of the low mass core simulation are shown in fig. II. H and 
fig. 11.21 they show the end state (at about 70,000 yrs) of the collapse of 
a B-E sphere that is chosen to be precisely the Bok globule observed by 



Alves et al.l (|200ll ) - whose mass is 2.1M© and radius R = 1.25 x 10 
AU at an initial temperature of 16 K. Two types of outflow can be seen: 
(i) an outflow that originates at scale of ~ 130 AU on the forming disk 
that consists of a wound up column of toroidal magnetic field whose 
pressure gradient pushes out a slow outflow; and (ii) a disk-wind that 
collimates into a jet on scale of 0.07 AU. A tight proto-binary system has 
formed in this simulation, whose masses are still very small < 1O -2 M0, 
which is much less than the mass of the disk at this time ~ 1O _1 M . 
Th e outer flow bears the ha llmark of a magnetic tower, f irst observed 
by lUchida k, Shibatal (|l985l ). and studied analytically by iLvnden-Bell 
d2003h . Both flow components are unbound, with the disk wind reaching 
3 km s" 1 at 0.4 AU which is quite super-Alfvenic and above the local 
escape speed. The outflow and jet speeds will increase as the central 
mass grows. 



In iBanerjee fc Pudritd (|2007h a similar setup to the one described 
above is used to study the formation of massive stars. In an unified pic- 
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ture of star formation one would also assume that outflows and jets will 
also be launched around young massive stars. Due to the deeply embed- 
ded nature of massive star formation observations of jets and outflows 
are much more difficult than for the more isolated low mass stars. Nev- 
ertheless, there is increasi ng observat i onal evidence for outflow s around 



young massive stars (e.g., Arce et al. . 2007 : Zhang et al. . 20071 ). 



The simulations of collapsing magnetised massive cores (~ 170 M Q ) 
show clear signs of outflows during the early stages of massive star for- 
mation. Fig. 11.31 shows two 2D snap-shots from this simulations where 
the onset of an outflow is visible. Early outflows of this kind might 
have a large impact on the accretion history of the young massive star. 



Krumholz et al J 420 051) showed that cavities blown by outflows help to 



release the radiation pressure from the newly born massive star which in 
turn relaxes radiation pressure limited accretion onto the central object. 

From these simulations one can conclude that the theory and compu- 
tation of jets and outflows is in excellent agreement with new observa- 
tions of many kinds. Disk winds are triggered during magnetised collapse 
and persist throughout the evolution of the disk. They efficiently tap ac- 
cretion power, transport a significant part portion of the disk's angular 
momentum, and can achieve different degrees of collimation depending 
on their mass loading. 



1.5 Feedback from collimated protostellar jets? 

The interstellar medium (ISM) and star forming molecular clouds are 



Elmegreen & Scalo, 


2004.1 Mac Low & Klessen, 


2004 


Ballesteros-Paredes et al., 


2007, and references therein). Supersonic turbulence 


in molecular clouds 



is known to play at least two important roles: it can provide pressure 
support to help support molecular clouds against rapid collapse, and 
it can also produce the system of shocks and compressions throughout 
such clouds which fragment the cloud into the dense cores that are the 
actual sites of gravitational collapse and star formation. One of the 
major debates in the literature is on the question of just how long this 
turbulence can be sustained - and with it, star formation. Despite the 
importance of supersonic turbulence for t he pr ocess of star formation 



its origin is still unclear. iNorman fc Silkl (j 19801 ) proposed that Herbig 



Haro (HH) outflows could provide the energetics to drive the turbulence 
in molecular clouds that could keep star formation active for a number 
of cloud free-fall times (i.e., several million years). This is an attrac- 
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tive idea as this proc ess could lead to a self-regulating star forma tion 
environment (see also Li fc Nakamural l2006t iNakamura fc Lil . l2007h . 

In a recent investigation by Ir3aneriee et al.l |2007t ) this idea was ad- 
dressed within a detailed study of feedback from collimated jets on their 
supra-core environment. These studies show that supersonic turbulence 
excited by collimated jets decays very quickly and does not spread far 
from the driving source, i.e. the jet. Supersonic fluctuations, unlike 
subsonic ones, are highly compressive. Therefore, the energy deposited 
by a local source, like a collimated jet, stays localised as the compressed 
gas either heats up (in the non-radiative case) or is radiated by cooling 
processes. The re-expansion of the compressed regions excites only sub- 
sonic or marginally supersonic fluctuations. This is contrary to almost 
incompressible subsonic fluctuations. These fluctuations travel like lin- 
ear waves with little damping into the ambient medium and make up 
most of the overall velocity excita ti ons. 

In the study of iBaneriee et al.l (|2007l ) a series of numerical experi- 
ments were performed with the adaptive mesh refinement (AMR) code 
FLASH. Here the jet is modelled as a kinetic energy injection from the 
box boundary. The energy injection could be switched on and off after 
a certain amount of time, i.e. the jets are either continuously driven 
or transient. These jets are either interacting with a homogeneous or 
clumpy environment, where the latter is modelled as a spherical over- 
density. Additionally the influence of magnetic fields on jet-excited fluc- 
tuations where considered in this study. 

The authors quantify the jet-excited motions of the gas using velocity 
probability distribution functions (PDF) which can be regarded as vol- 
ume weighted histograms. From these PDF the distinction between sub- 
sonic and supersonic fluctuations comes about naturally: the supersonic 
regime is strongly suppressed compared to the subsonic regime. Addi- 
tionally, separate decay laws of the kinetic energy for the subsonic and 
supersonic regimes where derived to quantify their different behaviours 
with time. 

Fig. 11.41 shows that the kinetic energy of the supersonic fluctuations 
decays much faster than subsonic contributions. Furthermore, one can 
see from fig. 11.51 that the jet-excited fluctuations are mainly subsonic 
which is due to the fact that the supersonic excitations do not travel far 
from the edge of the jet (see right panel of fig. 1 1 . 5 j) . 

This study shows that supersonic fluctuations are damped quickly 
because they excite mainly compressive modes. The re-expansion of 
the compressed overdensities drives mainly subsonic velocity fluctua- 
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Fig. 1.4. Time evolution ol the kinetic energies ,-Ekim f° r a transient jet with 
a speed ol Mach 5. The quantities are divided into a supersonic regime, v > c, 
and a subsonic regime, v < c. The decay ol supersonic energy contributions is 
much laster (laster than oc t~ 2 ) than the subsonic one. The time shift between 
the peaks in the kinetic energies shows that the subsonic fluctu ations are pow- 
ered b y the decay ol the supersonic motions. [Adapted from iBaneriee et al.l 
(2007)] 



tions that then propagate further into the ambient medium. This is in 
spite of the appearance of bow shocks and instabilities. In particular, 
instabilities such as Kelvin-Helmholtz modes at the edge of the jet, de- 
velop most efficiently for transonic or slower velocities. High-velocity 
jets, on the contrary, are bullet-like and stay very collimated, transit- 
ing the surrounding cloud without entraining much of its gas. From 
the point of view of jet-driven supersonic turbulence in molecular clouds 
this is a dilemma which is difficult to circumvent. Even in the case of 
overdense jets which affect more gas of the surrounding media and have 
higher momenta, the supersonic motions do not propagate far from their 
source. 

Simulations of magnetised jets in this study show that jets stay nat- 
urally more collimated if the magnetic field is aligned with the jet axis 
and therefore entrain less gas. Furthermore, perpendicular motions are 
damped by magnetic tension preventing a large spread of high ampli- 
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Fig. 1.5. Velocity structure and velocity PDFs at different times of a con- 
tinuously driven Mach 10 jet. Essentially no supersonic fluctuations in the 
ambient gas g et excited by the jet (the peak at v/c = 10 is the jet itself). 
[Adapted from iBaneriee et al.1 (120071 )] 



tude fluctuations. On the other hand, perpendicular field configurations 
support the propagation of such modes which are able to spread into a 
large volume. Nevertheless, the vast majority of these motions are still 
subsonic. 

This study shows that collimated jets from young stellar objects are 
unlikely drivers of large-scale supersonic turbulence in molecular clouds. 
Alternatively it could be powered by large scale fl ows which might be re- 



spons ible for the formation of the cloud itself (e.g.. lBallesteros-Paredes et al 



1999). Energy cascading down from the driving scale to the dissipation 



scale will then produce turbulent dens ity and velocity structure in the in- 
ertial range in between (|Lesieurl . ll997l ). If the large-scale dyna mics of the 



interstellar medium is driven by gravity (as suggested, e.g., bv lL i et al 



20051: iLi et all l200d) gravitational contraction would also determine to 
a large extent the internal velocity structure of the cloud. Otherwise, 
blast waves and expanding shells from super novae are also viable can- 
didates to power superson ic turbulence in molecular clouds (see e.g., 



Mac Low fc Klessenl . [2004) . 



1.6 Relativistic jets: theory 

The theory of relativistic MHD jets is best described in the SRMHD 
regim e. The essential parameter is the magnetisation parameter ([Michell . 
19691) . 



(1.18) 

4Mc 3 
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The iso-rotation parameter ilp is frequently interpreted as the angular 
velocity of the magnetic field lines. The function if? = B p r 2 is a measure 
of the magnetic field distribution (see Li 1993), and M = irpv p R 2 is the 
mass flow rate within the flux surface. Equation (|1.18[) demonstrates 
that the launch of a highly relativistic (i.e. highly magnetised) jet es- 
sentially requires at least one of three conditions - rapid rotation, strong 
magnetic field and/or a comparatively low mass load. 

In the case of a spherical outflow ("f? = const) with negligible gas pres- 
sure one may derive the Michel scal ing between the asymptotic Lorentz 
factor and the flow magnetisation (jMichell . Il969l ) , 



T - cr 1/3 



(1.19) 



Depending on the exact magnetic field distribution ^>(r,z) (which de- 
scribes the opening of the magnetic flux surfaces), in a collimating jet 
the matter ca n be substantially accelerated beyond the fast point magne- 
tosonic point (jBegelman fe Lil . ll994HFendt fc Camenzindl . ll996l ). as it is 
moved from infinity to a finite radius of several Alfven radii. As a result, 
the power law index in Eq . (11.191) can be different from the Miche l-scaling 
(see lFendt fc Camenzind Il996l ; Iviahakis fc Konigl bOQll 120031 ) . 

The light cylinder (hereafter i.e.) is located at the cylindrical radius 
ri = c/ftp. At the I.e. the velocity of the magnetic field lines "rotating" 
with angular velocity Qf coincides with the speed of light0. The I.e. 
has to be interpreted as the Alfven surface in the limit of vanishing 
matter density (force-free limit). The location of the I.e. determines 
the relativistic character of the magneto sphere. If the light cylinder is 
comparable to the dimensions of the object investigated, a relativistic 
treatment of MHD is required. 

While SRMHD can be invoked to study jets at scales larger than 
the gravitational radius of the black hole, close to the horizon general 
relativity becomes relevant and the MHD equations need to be coupled to 
General Relativity. This often requires the use of sophisticated numerical 
codes in order to capture the complexity of jet physics in GRMHD. In 
the (non-rotating) black hole's Schwarzschild spacetime the GRMHD 
equations are identical to the SRMHD equations in general coordinates, 
except for the gravitational force terms and the geometric factors of the 
lapse function. 



f Outside the I.e. the magnetic field lines "rotate" faster than the speed of light. 
As the field line is not a physical object, the laws of physics are not violated. 
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Fig. 1.6. Lorentz f actor from GRMHD simulations; adopted from 
iDe Villiers et a l, (2005). In all panels, the axes are in units of M (here 1M fa 4 
km). The black hole is located at the origin. The top row, left to right, shows 
plots of Lorentz factor for decreasing black hole rotation a/M = 0.995, 0.9, 0.0, 
respectively. The dotted contour marks the boundary of the jets. The only 
region where elevated values of T are found is in the jets (and also in the bound 
plunging inflow near the black hole, which is not resolved at the scale of this 
figure). Maximum values of V are found in knots that appear episodically in 
the upper and lower parts of the funnel. The bottom three panels show, from 
left to right, the corresponding time-averaged value of the Lorentz factor. The 
plots show evidence of spin-dependent collimation: cylindrical collimation is 
seen in the high-spin models, while the zero-spin model shows no such col- 
limation. These plots also show evidence of an extended acceleration zone: 
large Lorentz factors are built up over the full radial range. 



1.7 SRMHD and GRMHD Simulations 



Full GRMHD numerical simulati ons on the formation of jets near a black 
hole were first been performed by Kudoh et al. ( 1998h . This was followed 
by a plethora of GRMHD codes with fixe d spacetimes used to investi- 
gate jets from accreting black holes (e.g.. IDe Villiers fc Hawlevl. 12 003: 
Gammie et al. , 20031 Komissarov , 2004 : Anton et al. . l2006l : lAnninos et al. . 



2005). Here we focus on the most recent simulations that have managed 
to reach the highest Lorentz factors for the study of GRMHD jets. 
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Fig. 1.7. Asymptotic Lorentz factor (u v ,oo ) from SRMHD simul ations, for 
different magnetisation parameter a (from [Fendt fc Quvedl [2004). Lorentz 
values as high as 5000 can be reached for strongly magnetised sources. 



Recent state-of-the-art GRMHD simulations have shown that the ac- 
cretion flow launches energetic jets in the axial funnel region (i.e. low 
densit y region) of the disk /je t system, as well as a substantial coronal 



wind (|De Villiers et all I2005T ). The jets feature knot-like structures of 



extremely hot, ultra-relativistic gas; the gas in these knots begins at 
moderate velocities near the inner engine, and is accelerated to ultra- 
relativistic velocities (Lorentz factors of 50, and higher) by the Lorentz 
force in the axial funnel. The increase in jet velocity takes place in an 
acceleration zone extending to at least a few hundred gravitational radii 
from the inner engine. The overall energetics of the jets are strongly spin- 
dependent, with high-spin black holes producing the highest energy and 
mass fluxes. In addition, with high-spin black holes, the ultra-relativistic 
outflow is cylindrically collimated within a few hundred gravitational 
radii of the black hole, whereas in the zero-spin case the jet retains a 
constant opening angle of approximately 20 degrees. 

Figure 11.71 shows the Lorentz fact or, T, in the funnel ou tflow from 
the GRMHD simulations discussed in be Villiers et al.l|2005h . Elevated 
values of the Lorentz factor are found in compact, hot, evacuated knots 
that ascend the funnel radially; the combination of low density and high 
temperature in the knots^. The highest values of Lorentz factor reach 



f IQuved et al.l ||1997T ) discuss the appearance of knots in MHD jet simulations; the 
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the maximum allowed by the code (~ 50, and test runs show that much 
higher values could be reached with a higher ceiling) , and these are only 
found at large radii, suggesting the presence of an extensive region where 
the knots are gradually accelerated to higher Lorentz factors. 

While the details of the mechanism by which GRMHD jets are launched 
are yet to be fully understood, the following summary addresses the 
salient points of the complex dynamics by which un bound the plasma 



from the vicinity of the compact sta r (see discussion in lDe Villiers et al 



20051 : McKinnev fc Naravanl 2007 ): In the region of the accretion flow 



near the marginally stable orbit, both pressure gradients and the Lorentz 
force act to lift material away from the equatorial plane. Some of this 
material is lau nched magneto-cent r ifugal ly in a manner reminiscent of 
the scenario of iBlandford fc Pavnd (| 19820 . generating the coronal wind; 
some of this material, which has too much angular momentum to pen- 
etrate the centrifugal barrier, also becomes part of the massive funnel- 
wall jet. There is also evidence that the low-angular momentum fun- 
nel outflow originates deeper in the accretion flow. Some of this mate- 
rial is produced in a gravitohydromagnetic interaction in the ergosphere 
( Punslv fc Coronitil.ll990l). an d poss ibly in a process similar to that pro- 
posed bv lBlandford fc Znaiekl (jl977l ) where conditions in the ergosphere 
approach the force- free limit. The material in the funnel outflow is ac- 
celerated by a relatively strong, predominantly radial Lorentz force; gas 
pressure gradients in the funnel do not contribute significantly. 

In the context of GRBs, the Lorentz factor V of the relativistic wind 
must reach high values (P ~ 10 2 — 10 4 ) both to produce 7-rays and to 
avoid photon-photon annihilation along the l ine of sight, who se signa- 
ture is not observed in the spectra of GRBs ( Goodman . 198(jh . As we 
have seen, GRMHD simulations indicate Lorentz factors of up to ~ 50 
close to the black hole with indication of higher T far beyond the cen- 
tral object where SRMHD treatment is sufficient. The extreme Lorentz 
factors have been confirmed by recent SRMHD simulations which have 
fo und solutions with 7 as high as ~ 5000 as shown in Figure 11.71 (see 



Fendt fc Ouvedl |2004 for more discussion). 



knots seen in those simulations consisted of high density material, in contrast to 
what is observed in the present simulations. 
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1.8 Non-relativistic versus relativistic MHD jets 

While SRMHD and GRMHD jet simulations show similarities with non- 
relativistic regime, there are nevertheless some important differences in 
the underlying physics. 

• The existence of the I.e. as a natural length scale in relativistic MHD 
is not consistent with the assumption of a self- similar jet structure 
(as is often assumed in non-relativistic MHD jet models). The latter 
holds even more when general relativistic effects are considered. 

• Contrary to Newtonian MHD, in the relativistic case electric fields 
cannot be neglected. The poloidal electric field component is directed 
perpendicular to the magnetic flux surface. Its strength scales with 
the I.e. radius, E p = E± = (r/ri)B p . As a consequence of E p ~ B p , 
the effective magnetic p ressure can be lowered by a substantial amount 



(jBegelman fe Lil . 1 1994} 



In relativistic MHD the poloidal Alfven speed ma becomes complex 
for r > n, u\ ~ Bp (1 - (r/r ; ) 2 ) = b% — E\ . Therefore, Alfven waves 
cannot propagate beyond the I.e. and only fast magnetosonic waves 
are able to exchange information across the jet. 

In the relativistic case, break-up of the MHD approximation is an 
issue. The problem is hidden in the fact that one may find arbitrarily 
high velocities for an arbitrarily high flow magnetisation. However, 
an arbitrarily high magnetisation may be in conflict with the intrinsic 
MHD condition which requires a sufficient density of ch arged particle s 
in order to be able to drive the electric current sy stem (Michel . 19691) . 



Below the Goldreich- Julian particle density no (jGoldreich fe Julian! . 



19691 ) the concept of MHD as applied to relativistic jet breaks down. 



The points above imply that with due regard to the break down of 
the MHD approximation, it is feasible to scale the physics of jets from 
non-relativistic to ultra-relativistic MHD regime. GRMHD simulations 
suggest that a continuous scaling is more likely for slowly rotating black 
holes (where the radius of marginal stability is rather large and in a 
comparatively low-gravity region of spacetime) but show that a full gen- 
eral relativistic electrodynamic treatment is required to robustly treat 
the case of rapidly rotating black holes. GRMHD simulations also show 
that the situation is even more complex in the "plunge region" . This is 
the region of the disk within the radius of marginal stability in which 
the accretion flow is undergoing rapid inwards acceleration (ultimately 
crossing the event horizon at the velocity of light as seen by a locally 
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non-rotating observer). Unless the magnetic field is extremely strong, 
this is a region where inertial forces will dominate and the commonly em- 
ployed force-free approximation breaks down. As a particle gets closer 
to the central object it crosses three regions defining the MHD, force-free 
and inertial regimes, respectively. Another crucial aspect of the physics 
inherent to ultra-relativistic jets emanating from the near vicinity of 
rapidly rotating central obj ects is the role played by the dynamically im - 



Lay 

ed, 



portant electric field (e.g.. lLvutikov fc Ouvedl . 120071 : lOuved et all 120071 ) 
in driving plasma instabilities generally leading to pair creation. Pair- 
creation and subsequent annihilation into radiation has yet to be taken 
into account, at least consistently, in GRMHD or/and GREMD codes 
which are mostly based on mass conservation schemes. 

In summary, it appears that a universal aspect of all jets is their abil- 
ity to tap their energy from the accretion of gas into the gravitational 
potential of the underlying central object. Jets of all stripes may also 
share common morphological features far beyond the source. However, 
complete unification of non-relativistic with relativistic regimes for jets 
probably breaks down close to the central black holes, where jet physics 
will also depend on the rotation of the hole and the production of rela- 
tivistic plasmas. 
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